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ABSTRACT 

We study observed correlations between supermassive black hole (BHs) and the properties of their host 
galaxies, and show that the observations define a BH "fundamental plane" (BHFP), of the form M BH oc 
cr3-Oio.3_gO.43io.i9 or MjJH ^ _j_f0.54io.i7 cr 2.2±o.5 ) ana i ogous t0 the FP of elliptical galaxies. The BHFP is pre- 
ferred over a simple relation between Mbh and any of a, M*, Md yn , or R e alone at > 3 a (99.9%) significance. 
The existence of this BHFP has important implications for the formation of supermassive BHs and the masses 
of the very largest black holes, and immediately resolves several apparent conflicts between the BH masses 
expected and measured for outliers in both the Mbh - o and Mbh - M* relations. 
Subject headings: quasars: general — galaxies: active — galaxies: evolution — cosmology: theory 
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1. INTRODUCTION 

Discoveries of correlations between the masses of su- 
permassive black holes (BHs) in the centers of nearby 
galaxies and the properties o f their host spheroids (e.g., 
iKormendv & Richstond 1 19951) demonstrate a fundamental 
link between the growth of BHs and galaxy formation. A 
large number of similar correlations have now been identified, 
linkin g BH ma ss to host luminosity ([K ormendv & Richstond 
1 1995b. mass (jMagorrian et al.l |1998[). velocity dispersion 
dFerrarese & Merritj|2000t iGebhardt et alJl2000l). concentra- 
tion o r Sersic index (Graham et al.l l200lF iGraham & Driverl 
2006), and binding energy dAller & Richstond 120071) . among 
others. However, because these properties of host spheroids 
are themselves correlated, it is not clear whether any are in 
some sense more basic (see e.g. lNovak et al.ll2006l for such a 
comparison). 

The lack of a clear motivation for favoring one rela- 
tion over another has led to substantial observational and 
theoretical debate over the "proper" correlation for sys- 
tems which may not lie on the mean correlation between 
host properties, a nd over the demographics of the most 
massive BHs (e.g.. iBernardi et al.ll2006t lLauer et al.l l2006bt 
Batcheldor et al. 20061: IWvithdl200aT One possibility is that 
these different correlations are projections of the same "fun- 
damental plane" (FP) relating BH mass with two or more 
spheroid properties such as stellar mass, velocity dispersion, 
or effective radius, in analogy to the well-established funda- 
mental plane of spheroids. For the case of spheroids, it is 
now understood t hat various correlations , including the Faber- 
Jackson relation (Faber & Jackson 1 1976b between luminos- 
ity (o r effectively stella r mass M*) and velocity dispersion 
cr, the iKormendvl (1 19771) relation between effective radius R e 
and surface brightne ss I e , and the si ze-luminosity or size- 
mass relations (e.g., IShen et all 120031) between R e and M*, 
are all projections of a fundamental plane relatin g R e oc a a Ig 
dDressler et al.lll987l:lDlorgovski & Davislll987l) . 

In their analysis of the relation bet ween BH mass and hos t 
luminosity or dynamical mass, Mn y jMarconi & H unt (2003) 
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(see also Ide Francesco et all [2006) noted th at the residuals 
of the Mbh _ o relation (effectively Mbh/c 4 ; Tremaine et al. 
l2002l) were significantly correlated with the effective radii of 
the systems in their sample. Figure [TJ shows this, for the com- 
pilation of observations described in § [2]- there is indeed a 
clear trend that systems with larger effective radii or stellar 
masses tend to have larger Mbh/ (Mbh(c)). In that context, 
the authors argued for this as evidence favoring a relation be- 
tween M B h and Md yn oc cr 2 R e over Mbh oc ct 4 , but it is not 
clear that a dependence on Md yn alone completely or accu- 
rately captures the behavior in these residuals (and in §[3] we 
show that it does not). Furthermore, finding Mbh/o -4 oc R@ 
does not necessarily imply a FP-like relation, if the corre- 
lation between Mb h and a or a and R e has some nonlinear 
(e.g.,|\Vyithe 2006) or otherwise incompletely accounted-for 
behavior. Still, this brings up the important possibility of a 
true FP-like relation in which the combination of two proper- 
ties such as M* and a determines Mbh, which we study herein. 

Throughout, we adopt a £1m = 0.3, £1a = 0.7, Hq = 
70kms~ 1 Mpc~ 1 cosmology (and correct all observations ac- 
cordingly), but note this choice has little effect on our conclu- 



2. THE DATA 

We consider the sample of local BHs for which masses 
have been reliably determined via either kinematic or maser 
measurements. Specifically, we adopt the sample of 38 lo- 
cal systems for which values of Mbh, cr, Re, Mj V I | , and 
bulg e luminosities are co mpiled in Marconi & Hunt (120031) 
and lHaring & Rixl (J2004l) (see al so Mago rrian et al.l [19981 : 
iMerritt & Ferraresa 1200 U iTremaine et alj|2002l) . We adopt 
th e dynamical masses f rom the more detailed Jeans modeling 
in lHaring & Rix (2004). We estimate the tot al stellar mass M, 
from t he total /T-band luminosity given in iMarconi & Hunt! 
(120031 using the A _-band ma s s-to-li ght ratios as a function 
of luminosity from iBell et ail (120031) (specifically assuming 
a "diet" Salpeter IMF, although this only affects the abso- 
lute normalization of the relevant relations). We have re- 
peated these calculations using mass-to-light ratios estimated 
independently for each object from their B — J colors or 
full UBVJHK photometry, and find it makes little differ- 
ence. Where possible, we upda te measurements o f R e an d 
a with more recent v alues from Lauer et al. (2005, 2006a); 
iMcDermid et all (120061) and from lKormendy et all (120071) . 

Although it should only affect the normalization of the re- 
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FIG. 1 . — Observed residual in BH mass relative to the best-fit M^u — a re- 
lation, as a function of the host galaxy effective radius R e , velocity dispersion 
cr, and stellar mass M» . There is a clear suggestion in the data of a residual 
correlation with R e and M t , which we examine in detail in Figure|2] 

lations herein, we note that our adopted cosmology is identi- 
cal to that used to determine all quoted values in these works. 
When we fit the observations to e.g. the mean Mbh~c relation 
and other BH-h ost relations, we consid er only the subsample 
of 27 objects in lMarconi & Hunt! d2003) which are deemed to 
have 'secure' BH and bulge measurements (i.e. for which the 
BH sphere of influence is clearly resolved, the bulge profile 
can be well-measured, and maser spots (where used to mea- 
sure Mbh) are in Keplerian orbits). Our results are not quali- 
tatively changed if we consider the entire sample in these fits, 
but their statistical significance is somewhat reduced. 

3. A BLACK HOLE FUNDAMENTAL PLANE 

We wish to determine whether or not a simple one-to-one 
correlation between e.g. Mbh and a is a sufficient descrip- 
tion of the data, or if there is evidence for additional depen- 
dence on a second parameter such as R e or M*. The most 
efficient way to determine such a dependence is by looking 
for correlations between the residuals of the various projec- 
tion s of such a potent i al BH FP relation. In Figure Q] we fol- 
low iMarconi & Huntl (120031) and plot the dependence on the 



residual of the Mbh _ cr relation on R e and M* . This provides a 
clear suggestion of a residual dependence, and motivates us to 
examine a more complete description than a simple Mbh _ cr 
relation. However, in order to robustly do so, we must con- 
sider the correlations between residuals at fixed cr, not the cor- 
relation between e.g. the residual of Mbh _ cr and the actual 
value of a quantity such as R e or M* . 

Figure [2] plots the correlation between BH mass Mbh and 
host bulge effective radius R e or bulge stellar mass M*, all at 
fixed a. Specifically, we determine the residual with respect 
to the M B h - cr relation by fitting M B h(ct) to an arbitrary log- 
polynomial 



(log(M BH )) = S [a„ log(cr)"] 



(1) 



allowing as many terms as the data favor (i.e. until A% 2 with 
respect to the fitted relation is < 1), and then taking 

Alog(M BH | cr) ee log(M BH )- (log(M BH ))(<T). (2) 

We determine the residual A \og(R e | a) (or A log(M* | cr), for 
the stellar mass) in identical fashion, and plot the correla- 
tion between the two. We allow arbitrarily high terms in 
log(cr) to avoid introducing bias by assuming e.g. a simple 
power-law correlation between Mbh and cr, but find in prac- 
tice that such terms are not needed - as discussed below, 
there is no significant evidence for a log-quadratic (or higher- 
order) dependence of Mbh on a, R e , or M„, so allowing for 
these terms changes the residual best-fit solutions in Figure|2] 
by <C lcr. For reference, we find (in this sample) a best-fit 
correlation between R e oc a' i 0±0 - 4 (ignoring the log-quadratic 
and higher order terms, which are no t significan t ), sim ilar to 
that found for all SDSS galaxies by IShen et all (120031) . and 
M* oc cr 3 - 4±0 - 5 , again similar to the Faber- Jackson relation es- 



timated for global populations (e.g. lBernardi et al. 2003a|). In 
any case, changing these slopes slightly (or e .g. using the 
Mbh _ cr relation of Ferrarese & MerritJd2000h, as oppose d 
to our best-fit which is very similar to lTremaine et alj (2002)) 
has a relatively small systematic effect on the best-fit residual 
slopes (this generally makes the trend in residuals stronger, 
because the "mean" relation being subtracted is slightly less 
accurate for this particular sample), but does not qualitatively 
change our conclusions. 

Of course, even this approach could in principle introduce 
a bias via our assumption of some functional form, and so 
we have also considered a non-parametric approach where 
we take the mean (log(MBH)} in bins of log(cr). Unfortu- 
nately, the small number of observations limits such a non- 
parametric technique and somewhat smears out the interesting 
correlations (slightly decreasing their statistical significance). 
Regardless, the lack of evidence for a non-linear (i.e. log- 
quadratic or higher-order) relation means that we find similar 
results in all these cases, so we conclude that our methodol- 
ogy in determining residuals is not introducing a significant 
bias. 

The figure demonstrates that there is a highly significant 
correlation between Mbh and R e or M* at fixed a. We repeat 
this exercise in the figure, and demonstrate similarly that there 
is a highly significant correlation between Mbh and a or M* at 
fixed R e , and between Mbh and a or R e at fixed M*. This in- 
dicates that a simple one-variable correlation (e.g. a Mbh(ct), 
or Mbh(M„), Mbh(^c) relation) is an incomplete description 
of the observations. We therefore introduce a FP-like relation 
of the form 

M BH «cr Q flf, (3) 
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FIG. 2. — Left: Observed residual in BH mass as a function of host galaxy effective radius R e {top) or stellar mass M» (bottom), at fixed velocity dispersion 
a (equivalently, correlation between the residuals in the Mbh — a and R e -cr or M„ — a relations at each cr). The fit to this residual correlation is shown with the 
black lines (± 1 cr range in the best-fit correlation shown as dashed lines - note that they are strongly inconsistent with zero correlation), with the slope shown 
(dotted blue line shows the least-squares bisector). The probability of the null hypothesis of no correlation in the residuals (i.e. no systematic dependence of Mbh 
on R e or M* at fixed cr) for the observed systems is shown (red P nl ,n). Center: Residual in Mbh as a function of cr or R e at fixed stellar mass M* . Right: Residual 
in Mbh as a function of cr or M* at fixed effective radius R e . The observations imply a secondary "fundamental plane"-type correlation at 3cr with respect to 
each of these variables. 



which can account for these dependencies. Formally, we de- 
termine the combination of (a, (3) which simultaneously min- 
imizes the x 2 / ly of the fit and the significance of the cor- 
relations between the residuals in rj and Mbh (or R e and 
M bh)- This yields simila r results to the direct fitting method 
of iBernardi et alJ (12003b) from the spheroid FP, which mini- 
mizes 



TABLE 1 
BH-Host Correlations 



A 2 = [ log(M BH ) - a log(a) - /3 logCR.) - S] ! 



(4) 



It is straightforward to extend this minimization by weighting 
each point by the measurement errors (where we allow for 
the errors in all observed quantities - log(MBH), log(cr), and 
log(/? e ), and estimate symmetric errors as the mean of quoted 
two-sided errors). This yields a best-fit BHFP relation 



Variables ' 


Normalization 2 


Slope 3 


Scatter 4 


a a R e 3 


8.33 ±0.06 


3.00 ±0.30, 0.43 ±0.19 


0.21 




8.24 ±0.06 


0.54±0.17, 2.18±0.58 


0.22 


M?fif 


8.06 ±0.07 


1.78 ±0.40, -1.05 ±0.37 


0.25 


M,a 2 


8.23 ±0.06 


0.71 ±0.06 


0.25 


a 


8.28 ±0.08 


3.96 ±0.39 


0.31 


M» 


8.21 ±0.07 


0.98±0.10 


0.33 


M dyI1 


8.22 ±0.10 


1.05±0.13 


0.43 


Re 


8.44 ±0.10 


1.33±0.25 


0.45 



log(M BH ) = 8.33 + 3.00(±0.30) log(cr/200km S - 1 ) 
+0.43(±0. 19) log(fl e /5 kpc) 



(5) 



from the observations. Unsurprisingly, the slopes in the BHFP 
relation are close to those formally determined for the residu- 
als in Figure|2] As expected, the residuals of Mbh with respect 
to these fundamental plane relations, at fixed R e and fixed 
a, show no systematic trends and are consistent with small 
intrinsic scatter. The introduction of a BHFP eliminates the 
strong systematic correlations between the residuals, yielding 
flat errors as a function of a and R e . 

At low redshift, cr, R e , and Md yn can be determined reli- 
ably, but at high redshift it is typically the stellar mass M* or 
luminosity which is used to estimate Mbh- Therefore, it is 
interesting to examine the BHFP projections in terms of e.g. 
M* and a or M* and R e . Repeating our analysis, we find in 
Figure [2] that the observations demand a FP relation over a 
simple Mbh(M») relation at high significance. The exact val- 
ues of the best-fit coefficients of this BHFP determined from 



For the variables (x, y), a correlation of the form log(MBH) = ct log(x) + log(y) + 5 
is assumed, where the normalization is S and a, {3 are the logarithmic slopes. The 
normalization gives log(M BH ) for cr = 200kms~', M, = lO n M , M dyn = 10 1! M o , 
R e = 5 kpc, which roughly minimizes the covariance between fit parameters. 3 Errors 
quoted here for the BHFP relations in {a, R e ), (M* . a), and (M* , R e ) include the covari- 
ance between the two slopes. Holding one of the two fixed and varying the other yields 
substantially smaller errors (typically ~ 5%). All quoted errors account for measure- 
ment errors in both Mbh an d the relevant independent variables. 4 The internal scatter is 
estimated from the observations as that which yields a reduced 'X 1 jv = 1 with respect to 
the given best-fit relation. 

the observations are given (along with those of various other 
BHFP projections) in Table Q] 

The BHFP in terms of a and R e (i.e. M BH oc a a R^) is 
of course tightly related to the BHFP in terms of a and M* 

(Mbh oc a a M* ) or M* and R e ; the near-IR fundamental plane 
relates stellar mass (assuming /T-band luminosity is a good 
proxy for stellar mass), effective radius, and velocity disper- 
sion as R e oc cr 1 53 If 19 (iPahreet al.ll998l) . Using/, cxM*/R 2 e , 
we can substitute this in Equation © and obtain the expected 
BHFP in terms of a andM*, namely M B h oc rj L9 M"- 58 . This is 
quite close to the result of our direct fitting, Mbh oc cr 2 - 2 M"- 54 . 
The FP of early-type galaxies, in other words, allows us to 
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relate our two "fundamental planes," namely Mbh oc a R e 

1 /2 

and Mbh oc a MJ . Given the tight early-type FP relation 
between M„, <j, and R e , these two forms of the BHFP are 
completely equivalent (the choice between them is purely a 
matter of convenience). We can, in fact, pick any two of 
the three FP-related variables as our independent variables for 
predicting Mbh - using the early-type FP again to transform 
the BHFP relation in terms of a and R e to one in terms of 
M» and R e , we expect M B h oc Ml 6 R~ S , similar to the rela- 
tion we directly fit (see Table [TJ. Any two of M*, R e , and a 
can thus be used to predict Mbh according to the BHFP rela- 
tions. The tightness of the early-type FP also means that it is 
redundant to search for a four-variable correlation (i.e. one of 
the form Mbh oc a a R^ M*), since M* is itself a function of 
a and R e with very little scatter (and indeed, directly testing 
this, we find no significant improvement in our fits expand- 
ing to a four-variable correlation). Perhaps most important, 
however, is that no transformation (given the early-type FP 
relating these three variables) eliminates the dependence on 
two variables - i.e. no transformation from any one BHFP al- 
lows us to write Mbh as a pure function of either a, R e , M*, 
or Mdyn (as we expect, since Figure [2] explicitly shows that 
each of these single-variable correlations exhibits a significant 
residual dependence on another variable). 

Given the definition of Md yn oc <r 2 R e , it is trivial to convert 
the best-fit BHFP relation in terms of a and R e to one in Md yn , 
obtaining 

log(M BH ) = 8 .29 + 0.43 log(M dyn / 10 1 1 M ) (6) 
+2.141og(cr/200kms- 1 ) 
= 8.19+ 1.50 log(M dyn /10 n M Q ) 
-1.07 log(fi e /5kpc). 

It is important to note that the residual correlation with either 
a or R e at fixed Md yn is non-zero and highly significant. It is 
therefore not the case that the BHFP reflects, for example, the 
true correlation being between Mbh and M<j yn (in which case 
the BHFP would have a form M BH oc (a 2 R e )t), but it is a FP 
in a genuine sense. 

We have also repeated the analysis of Figure |2]forMd y n and 
a (or Ma yn and R e ), and obtain these results directly, with a 
pure correlation between Mbh and Md yn ruled out at ~ 3a in 
the observations. We note that any analysis of these particular 
residual correlations must be done with special care, as the 
plotted quantities (e.g. MBH/Mi y n versus R e / ' {R e [M d yn ])) are 
not independent (since M dyn oc a 2 R e /G depends directly on 
the measured R e ), and so fitting a small sample where there is 
both intrinsic scatter and measurement errors in the quantities 
can bias the fit. Running a series of Monte Carlo experiments 
allowing for the range of estimated intrinsic scatter and mea- 
surement errors in each quantity, we find that the observations 
span a sufficiently small baseline that, with the present errors, 
a naive comparison will be biased to underestimate the sig- 
nificance of the preference for a BHFP relation over a simple 
correlation with Md yn . The P nu u for correlations between the 
residuals in the MBH~Md yn relation and R e or a typically de- 
creases by a factor ~ 2 with a proper Monte Carlo analysis; 
i.e. the true significance of the BHFP relation is even greater 
than a direct comparison suggests. 

Figure |2] demonstrates the significance with which the ob- 
servations rule out both a pure BH-host mass relation (either 
Mbh ocM» or Mbh ocMd yn ) and a pure Mbh -cr relation. How- 
ever, when fitting to a form Mbh oc a a R@, there is still some 




a = 4-4(3 (M BH oc CT a lvlP) 

FIG. 3. — Probability that there is no remaining trend or correlation among 
the residuals of the Mgn-host relation (P nl iii, as in Figure f2), as a func- 
tion of the slope a for a correlation of the form Mbh oc a a Rg {top) or 

Mbh oc o~ a Mf (bottom). For each value of a, we marginalize over j3 and the 
correlation normalization to find the best fit, but there is a rough degeneracy 
between the best fit a and j3 ((3 (4 - a)/2 for M B h oc o a R% , (3 ss (4 - a)/ 4 
for Mbh oc o a M* ). Lines show P nu n for the observations, points show the 
best fit (in a x 2 sense) BHFP relation and errors from Table [T] Thick black 
lines show the a corresponding to a pure Mbh ~o, Mbh ~^dyn, or Mbh — M* 
relation, all of which are ruled out at similar (~ 3<r) significance. 

degeneracy between the slopes a and (3 (roughly along the 
axis « (4 — a)/2). Figure [3] illustrates the degree of this 
degeneracy and the extent to which, for example, a BHFP 

"A 1/2 

with Mbh oc a R e is favored over a pure MBH-M dyn rela- 
tion. We plot the likelihood of a residual correlation between 
Mbh and R e or a at fixed a a R@, as a function of the slope a 
(marginalizing over (3 and other fit parameters at each a, al- 
though (3(a) roughly follows the axis of degeneracy above). 
In detail, this is the likelihood of a correlation between the 
residuals of the M B h oc a a R& relation and the residuals of the 
a(a a RP) and R e {cr a R^) relations, identical to the procedure 
used to determine the quoted P nu u in Figure [2] We compare 
the quoted best-fit slopes from Table Q] We then repeat this 
exercise for the alternate representation of the fundamental 
plane, M B h oc a a M*. 
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The analysis shown in Figure[3]agrees reasonably well with 
the x 2 best-fit expectations, and illustrate an important point: 
a pure relation between BH and host mass (either dynamical 
mass Mdyn or stellar mass M*) is ruled out at a significance 
level comparable to that with which a pure Mbh - cr relation 
is ruled out (~ 3 cr). Indeed, the preferred BHFP parameters 
are centered almost exactly mid-way between these two pre- 
viously proposed relations. 

However, there are possible correlations which cannot be 
clearly discriminated by the present data. Both a pure relation 
between BH mass and spheroid binding energy, of the form 
M BH °c (M» cr 2 ) 2/3 , as studied in lAller & R ichstone (2003), for 

1/2 9 

example, and a mixe d relation of the form Mbh oc MJ a (see 
Hopkins et al. 2007), are within the ~ ler allowed range of the 
data. It is worth noting that simply expanding the number of 
observed sources will not necessarily break these degenera- 
cies. Rather, to increase the constraining power of the obser- 
vations, a larger baseline is needed, including (in particular) 
a larger sample of objects which lie off the mean R e -a or 
M, - a relations (and thus extend the baseline in the residual- 
residual space which properly constrains the FP slopes). 

Finally, we noted above that Figure [2] is essentially un- 
changed if we consider residuals with respect to just linear 
(i.e. pure power law) MBH~cr-type relations, as a consequence 
of there being no significant evidence in our data for a log- 
quadratic or higher-order correlation. Allowing log-quadratic 
terms in our fundamental plane fits, we find a best fit to the 
observations of the form 

log(M BH ) = 8 .06 +(2.8 ± 0.4) a+ (0.48 ± 0. 1 8) f (7) 
-(2.1 ±2.3)ct 2 + (0.31 ±0.25) r 2 

(where a = log [cr/200kms _1 ] and r = log [/? c /3kpc]). The 
linear BHFP coefficients in a and R e are similar to those in 
Table Q] and their significance is not much changed - this il- 
lustrates that the FP behavior we find cannot simply trade off 
with or be equally well-represented by a log-quadratic depen- 
dence (i.e. one cannot eliminate the residual dependence of 
Mbh on R e at fixed a by adding a log-quadratic or higher- 
order term in cr). The log-quadratic terms are at most signifi- 
cant at the ~ ler level. This is also true if we add just one of 
the two log-quadratic terms - adding a log-quadratic term in 
just R e yields a coefficient (0.16 ±0.25) r 2 , and adding one in 
just a gives a coefficient (- 0.28 ± 2. 15) a 2 . This is similar to 
the finding of l Wyithel (120061) . who estimates < ler significance 
for the addition of log-quadratic terms in any of cr, Md yn , or 
M*. 

4. CONCLUSIONS 

We study the correlation between observed central BH mass 
and host galaxy properties, and find that the systems lie on a 
BH "fundamental plane" (BHFP), of the formM BH oc a 30 R° e 5 
or Mbh oc M"' 5-0,7 cr 1 - 5-20 , analogous to the FP of spheroids. 
Specifically, there are significant (at > 99.9% confidence) 
trends in the residuals of the Mbh _ cr relation with M* and 
R e at fixed cr, and likewise in the Mbh~M» relation (with a 
or R e at fixed M»). This provides a new paradigm for under- 
standing the traditional relations between BH mass and either 



bulge velocity dispersion or mass. These correlations (as well 
as those with other bulge properties such as effective radius, 
central potential, dynamical mass, concentration, Sersic in- 
dex, and bulge binding energy) are all projections of the same 
fundamental plane relation. Just as the Faber-Jackson relation 
between e.g. stellar mass or luminosity and velocity disper- 
sion (M* - cr) is understood as a projection of the more funda- 
mental relation between M„, cr, and R e , so too is the Mbh _ cr 
relation (Mbh oc cr 4 ) a projection of the more fundamental re- 
lation Mbh oc cr 3 R~ ,s . Recognizing this resolves the nature of 
several apparent outliers in the Mbh _ cr relation, which sim- 
ply have unusual a values for their stellar masses or effective 
radii, and eliminates the strong correlations between residu- 
als. 

Improved measurements of the host properties of sys- 
tems with well-measured BHs can significantly improve con- 
straints on the BHFP. As noted in Table Q] the present ob- 
servations demand a correlation of the form Mbh oc cr a M* 
over a simple correlation with either a or M* at > 3 cr con- 
fidence. Already, this puts strong constraints on theoretical 
models of BH growth and evolution - BH mass does not sim- 
ply scale with the star formation (stellar mass) or virial ve- 
locity of the host galaxy. However, there is still a substan- 
tial degeneracy between the slopes a and f3 (roughly along 
the axis (3 ~ 1 - a/4). For example, the existing data do not 
allow us to significantly distinguish a pure correlation with 
spheroid binding e nergy Mbh oc (M*cr 2 ) 2/3 , as detailed in 
lAller & Richstond (|2007) from the marginally favored rela- 
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tion oc MJ cr . Both suggest that the ability of BHs to self- 
regulate their growth must be sensitive to the potential well 
at the center of the galaxy (and therefore to galactic struc- 
ture), but the difference could reveal variations in the means 
by which BH feedback couples to the gas on these scales. 

Increasing the observed sample sizes and, in particular, ex- 
tending the observed baselines in mass and a will substan- 
tially improve the lever arm on these correlations. In partic- 
ular, the addition of stellar mass M* information to the sig- 
nificant number of objects which have measurements of cr 
and indirect measurements of Mbh from reverberation map- 
ping would enable considerably stronger tests of our proposed 
BHFP relation. We do note the caveat, however, that care 
should still be taken to consider only bulge properties and re- 
move e.g. rotationally supported contributions to the velocity 
dispersion. 

The BHFP appears to be a robust correlation, which pro- 
vides an improved context in which to understand the nature 
and evolution of the numerous observed correlations between 
BH and host spheroid properties. In particular, the results de- 
scribed here provide new, important constraints for models of 
BH growth, feedback and self -regulation. 
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that greatly improved the clarity of the manuscript. This 
work was supported in part by NSF grant AST 03-07690, and 
NASA ATP grants NAG5-12140, NAG5-13292, and NAG5- 
13381. 



APPENDIX 

DETERMINING THE SIGNIFICANCE OF RESIDUALS IN THE M BH -HOST CORRELATIONS 

We test the robustness of our results, and how appropriate our fitting method is, by constructing a series of Monte Carlo 
realizations of the observations. As the velocity distribution a is the most well-measured quantity for the observed galaxies (with 
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FIG. 4. — Top: Recovered dependence of Mbh on Re at fixed a (blue lines) from Monte Carlo realizations of the observed data sets, assuming different "true" 
intrinsic correlations between Mbh and host properties (of the form Mbh c< cr" Re)- Black lines in each panel show the residual of the assumed "true" relation, 

i.e. zero dependence in the case where the true relation is Mbh <x cr 4 (left), residual Mbh oc R l J 2 for a BHFP relation of the form Mbh oc a 3 R l J 2 (center), and 
residual Mbh oc Re for a Mbh oc Md yn relation (right). Solid blue lines show the median dependence recovered using our standard fitting method (fitting residuals 
at fixed a) for all Monte Carlo realizations, dashed lines show the 25/75% range of fits recovered. If we fit the residuals as a function of the actual value of 
Re (instead of R e at fixed cr), we infer the correlations shown as red lines. In all cases, fitting residuals versus residuals (our standard methodology) recovers a 
similar slope to the "true" correlation and appropriate (in a x 2 sense for normal errors) range about that true slope. However, comparing to the actual value of R e 
severely smears out the true dependence on R e when it is present. Bottom: The distribution of the significance of the residuals (of Mbh with respect to R e ) in our 
Monte Carlo realizations. The same trends are seen: fitting residuals versus residuals yields the expected \ 2 distributions and generally detects the appropriate 
significance of residuals when there is a true dependence on R e , but in almost all cases (even when the true correlation has a strong dependence on R e ), fitting 
residuals in Mbh against the actual value of R e leads one to conclude (incorrectly) that the residuals are not significant. 

the exception of the Milky Way, which we generally exclude from our fits following Tremai ne et al.l (120021) ). we begin with the 
observed distribution of a values. Either from the cumulative distribution of a values, or from the observed error bars for each 
value of cr, we then statistically resample the a distribution (with the same number of objects) for each Monte Carlo realization. 
Next, we assume a mean intrinsic correlation between R e and a (for simplicity, assume R e oc cr 2 , similar to that observed), with 
some intrinsic scatter comparable to that observationally inferred, and use this to randomly generate the true R e values of each 
point. Then, given some assumed true BH-host correlation of the form Mbh oc a a R^ , again with an intrinsic scatter comparable 
to the observational estimates, we randomly generate the true BH masses. Finally, using the mean observational measurement 
errors in each of these quantities (or, specifically the quoted observational errors of each object corresponding to each mock point, 
as it makes little difference), we randomly generate the "observed" values of each quantity. We then repeat our fitting procedures 
from § [3] For example, to search for residuals with respect to the M BH - a relation, we fit the "observed" points to a mean 
M BH - cr relation, and compare the residuals in BH mass to the residuals in R e of the "observed" R e - a relation. We determine 
the statistical significance of these residuals, and fit to determine the residual dependence (i.e. best fit "observed" Mbh oc a a R% 
relation). 

Figure E] shows the results of this analysis. Specifically, we show the median correlations between the residuals in Mbh and 
in R e , from ~ 1000 Monte Carlo realizations, and the 25/75% quartile ranges of the fitted correlations. We show this for three 
assumed cases - one in which the "true" correlation is a pure M B h - cr relation (M BH oc cr 4 ), one in which the true correlation 

i 1 12 

is similar to our best- fit BHFP relation (Mbh oc a R e ), and one in which the true correlation is a pure MBH-Md yn relation 
(Mbh oc Md yn oc a 2 R e ). In all of these cases, our fitting method generally recovers a residual correlation very similar to the input 
correlation. There appears to be a (very slight) bias towards our fitting method recovering a marginally stronger dependence of 
Mbh on R e , however, in terms of its statistical significance, the bias introduced is very small - at most, it implies that the "true" 
(in a maximum likelihood sense) significance of the dependence of M B h on R e at fixed a which we recover from the observations 
should be slightly lower than quoted (equivalent to increasing our error bars by ~ 10%). In other words, the observational data 
span a sufficient baseline, and the fitting method we use is sufficiently robust, that any "true" correlation of the form Mbh oc a a R@ 
should be recovered in the majority of cases. 
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We can also check the significance of the residuals (of Mbh with respect to R e ) from each of these Monte Carlo realizations, as 
shown in Figure@] In the case of a pure "true" Mbh - cr relation, the residuals should not be significant in most cases, and indeed 
we find that they are not. Their behavior is essentially exactly what is expected for a \ 2 distribution - i.e. in 5% of all cases, a 
95% significance level is assigned to the residuals, in 1% of all cases, a 99% significance level is assigned. The significance as we 
have determined it therefore carries its appropriate (expected) meaning and weight, and our general \ 2 analysis and assumption 

t 1 /2 

of normal errors in log-space is probably not misleading. In the case where the "true" correlation is of the form Mbh oc a R e or 
Mbh oc a 1 R e , the likelihood that the residuals will be significant increases dramatically (as it should), although we caution that 
there is still some (fairly large) probability that the residuals of a given Monte Carlo realization will not be especially significant 
(< 2cr significance). From comparison with the residuals with respect to a pure Mbh-c relation, this implies that the parameter 
space spanned by current observations is such that a significant detection of residuals does imply a significant probability of a 
true BHFP relation, but that a weak or non-detection of such residuals does not yet rule out such a relation. 

We caution, however, that this is only the case when comparing residuals against residuals. If we were, as in Figure Q] to 
simply consider the residual in Mbh - c as a function of the actual value of R e (instead of comparing Mbh and R e both at fixed 
c), we smear out the significance of any real residuals, as we would expect. The slope recovered (i.e. the inferred dependence 
of Mbh on R e ) is severely biased towards being too shallow for any non-zero dependence on R e , and in only ~ 1 % of cases will 
such a method recover a slope similar to the "true" intrinsic correlation. Looking at the significance of the residuals in this space, 
it is clear that this projection biases against detecting any significant residual dependence on R e . Even if the true relation were 
Mbh oc a 2 R e , looking at the significance of the residuals as a function of the actual value of R e would lead one to conclude in 
over 80% of cases that they were insignificant (< 2a formal significance). A proper analysis, however, would recover the true 
significance of the residuals in > 95-99% of cases. 

Finally, by repeating this analysis over the entire parameter space of "true" correlations (for example forms Mbh oc a a R$), and 
comparing with the data, we can use a similar Monte Carlo approach to determine a maximum-likelihood, best-fit BHFP relation. 
We have done so, and find nearly identical results to our standard fits, again indicating that the conclusions herein are robust to 
the exact likelihood calculation so long as the residuals are properly analyzed. Similarly, we can repeat this entire analysis for 
true correlations of the form Mbh oc a a Mi or Mbh oc M" R@, or in order to test the significance of residuals at fixed R e or fixed 
M*. In all cases, we recover very similar results, reinforcing the robustness of our conclusions. 
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